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We derive quasiclassical expressions for the three-body decay width and define the "preexponen- 
tial" coefficients for them. The derivation is based on the integral formulae for the three-body width 
obtained in the semianalytical approach with simplified three-body Hamiltonian [L.V. Grigorenko 
and M.V. Zhukov, arXiv:0704.0920vl]. The model is applied to the decays of the first excited 3/2~ 
state of 17 Ne and 3/2" ground state of 45 Fe. Various qualitative aspects of the model and relations 
with the other simplified approaches to the three-body decays are discussed. 
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I. INTRODUCTION 

The "true" two-proton decay [l] is the decay mode 
which is expected to be an ordinary phenomenon in the 
vicinity of the proton dripline [2| . This mode corresponds 
to a specific situation when one-proton emission is en- 
ergetically (due to the proton separation energy in the 
daughter system) prohibited and only the simultaneous 
emission of two protons is possible. The energy condi- 
tions are illustrated in Fig. [TJ more detailed discussion 
of the three-body decay modes could be found in 0, Q • 
From theoretical point of view this situation is a subset 
of the three-body Coulomb problem in the continuum. 
A consistent quantum mechanical three-cluster model of 
the phenomenon was developed in Refs. 0, H, 0] and 
applied to a range of the nuclear systems from 6 Be to 
66 Kr in papers [I @, 0, @] . In the works [1, E3] possible 
importance of the "true" two-proton decay phenomenon 
was demonstrated for astrophysical applications (in the 
form of the reverse two-proton radiative capture process) . 
Having in mind these applications, which could require 
precise result in certain cases, the semianalytical model 
with simplified three-body Hamiltonian was developed in 
our recent work [llj. The model allowed precise calibra- 
tions of three-body calculations for decay widths. 

On the other hand the semianalytical model [ll| can be 
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FIG. 1: Energy conditions for different modes of two-nucleon 
emission (three-body decay): true three-body decay (a), se- 
quential decay (b) . 



used for derivation of the quasiclassical formulae for the 
three-body decay width. Formulae of this class can be 
found in early papers on two-proton radioactivity [J, . 
They were used for qualitative estimates in our works 
[3, B S H) 0j E3| • We characterize them as quasiclas- 
sical as they are based on a certain factorization of the 
decay amplitude, which, in reality, is equivalent to exis- 
tence of classical trajectories for system propagation in 
the process of the decay. A quasiclassical model for the 
width was introduced in the papers [lj| for 12 O and 
then utilized in the series of works by Barker and Brown 
[H [H, [H, [H HI HI in somewhat modified form. It 
was characterized as "R-matrix approach" , due to for- 
mal similarity with two-body R-matrix formalism. 

The name of "R-matrix approach" can be, in certain 
sense, misleading as a derivation procedure is actually not 
based on the construction of the R-matrix on the "nuclear 
surface" . The latter has a very complicated shape in the 
three-body system due to pairwise final state interactions 
(FSI). Fig. [2] shows schematically the surface limiting the 
region of classically allowed motion in the nuclear interior 
space for pairwise distances between clusters r*y . Close 
to the origin this surface is approximated by the surface 
of constant hyperradius (ellipsoidal in this space). In 
the regions of final state interactions ry <C , rjy it is 



const. Possible values of radii due to trian- 
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gle conditions |r.;fc — r/y| < < + r^j are limited 
by solid angle represented by gray tetrahedron in Fig. 
[2] When particles penetrate through the Coulomb bar- 
rier they propagate predominantly along the "tunnels" 
of classically allowed regions (though the motion entirely 
in the classically allowed region is not possible due to the 
energy condition of Fig. [T]). 

In this work we derive the quasiclassical expression 
for the three-body decay width in the case of the ex- 
istence of narrow states in the subsystems and define the 
"preexponcntial" coefficients for the expression. The ob- 
tained expression has an important advantage compared 
to those used previously. We also make a comprehen- 
sive overview of quasiclassical three-body formulae used 





The function ipi on the left from R is obtained as a so- 
lution of the Schrodinger equation [normalized only by 
asymptotic condition 'F; = G{\. To the right of the 
matching point it can be written explicitly in the form 
correctly normalized both asymptotically (to G/) and in 
the internal region (to unity): 
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FIG. 2: Qualitative plots of the surface limiting the region 
of classically allowed motion in the nuclear interior for three- 
body problem, (a) One FSI and (b) three FSIs. Coordinates 
Tij are pairwise distances for three constituents. The gray 
tetrahedron limits the possible values of r;j . 



in different works and provide a qualitative analysis of 
them. 

The unit system % = c = 1 is used in the article. 



II. TWO-BODY CASE 

The inte gra l formula for the decay width in two-body 
case is uU [2311 



2L = ± 



dr<pi(k r r) (V-V)Mk r ,r) 



(1) 



where the WF %pi{k ri r) is the "quasibound" solution for 
Hamiltonian 

(H - E r )fa = (T + V - Er)fa = 

at the resonant energy E r . "Quasibound" means that the 
solution is matching the irregular (at the origin) Coulomb 
function Gi in the subbarrier region and it is normalized 
to unity in the finite region of radius R. The function 
tpi (kr) is the continuum WF of the auxiliary Hamiltonian 
H 

(H - E) Vl = (T + V — E)ipt = 



Neglecting the terms of the order Si we obtain 
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For small widths the relation Gi 
penetrability can be identified as 
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Then we obtain for the width in Eq. |T]) the standard 
R-matrix expression 



r = 2 1 2 P l {k r ,R) = 2j^ L 9 2 Pi(k r ,R) , 



(3) 



where the Wigner limit for the reduced width r )^ VL and 
the dimensionless reduced width (DRW) 6 2 are identified 
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The DRW is expected to be of the order of unity for R 
chosen in the subbarrier region. However, the DRW 9 2 
decreases relatively strongly with radius R for two rea- 
sons: (i) the normalization integral in the denominator 
is growing with R [27J; (ii) the numerator of Eq. ((4| is 
decreasing as Gf with R under the barrier. This decrease 
is compensated in Eq.[3]by growth of the penetrability as 
1/Gf ■ Thus for values of R under the barrier the width 
should decrease slowly with radius R up to some value 
of the radius (close to the outer classical turning point of 
the barrier), where it is practically stabilized. 



in the S-matrix representation having the asymptotic 
form 

(fi(kr) = exp(£(Sj) [Fj(At) cos(<J;) + G t (kr) sin(5/)]l . 

_(2) 

The only restriction on the auxiliary Hamiltonian H is 
that its eigenfunction should be sufficiently far from res- 
onance at E = E r . 

The integral in E q. (U ) can be rewritten in terms of 
Wronskian, see Ref. [Ill ]: 
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A. Integration region in calculations of widths 

The role of the integration limit in the calculations of 
the internal normalization can be understood using the 
R-matrix example. Similar calculations can be found in 
the book 24j , where they illustrate a somewhat different 
issue. In the case of a square well potential and absence 
of Coulomb interaction the normalization integral for the 
internal region is easily evaluated 



[{kRx)ji{kRx)\ 2 dx = ^^[ji(kR) 2 



jt+tikRV^hR)] « MKR)] 2 ■ (5) 
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The approximate equality is good at resonance and be- 
comes exact if E r — >■ 0. The dimensionless reduced width 
is in that case 
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We can redefine the DRW in such a way that penetrabil- 
ity could be evaluated at the nuclear surface but the effect 
of the contribution to normalization from the subbarrier 
region is already taken into account. The normalization 
integral is reasonable to calculate up to the first zero 
of the irregular function Gi(k r r). The WF comes from 
under the barrier approximately at this point and cal- 
culation of the normalization integral beyond this point 
loose a sense (the particle can not anymore be considered 
to be in the "internal region"). 

Beyond the matching point R the function can be ap- 
proximated as 

Gi{k r r) = (k r r)ni(k r r) — r~ l , 

with a good precision and integration in Eq. (0)) can be 
extended to infinity. The integral of such approximated 
function converges for integration up to infinity for / > 0. 
Thus the normalization for the whole subbarrier domain 
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The "redefined" dimensionless reduced width is exactly 
written as 



B. Evaluation of the integral in Eq. (TTJ) 

If the width of the resonance in the auxiliary Hamilto- 
nian is sufficiently small, then, in proximity of the auxil- 
iary resonance energy E a , we can write confidently that 
in the internal region 



<Pi(kr) = \l-2\l 1 {E -E a y + Tl{E)/A 



i>i(k a ,r) , 



where function ^i(k a ,r) is also the quasibound WF as 
ipl(k r ,r) (namely with resonant boundary condition and 
normalized to unity in the internal domain) but taken at 
the resonant energy E a of the auxiliary Hamiltonian 



ipi(k a ,r > R) tx Gi(k a r) 



dr 



i>t(k a ,r) 



The coefficient Di provides the normalization of the 
Breit-Wigner profile, which is already practically normal- 
ized 



dE- 



T a {E a )/2ir 



(E-E a y + Tl{E a )/A 
within the energy domain of interest: 

^ +AE T a {E)/2-K 



= 1, 



l = Di 



dE 



{E - E a ) 2 + T a (E) 2 /4 



For reasonably narrow states Di = 1 means that AE = 
(3— 5)T a (E). So this coefficient is always reasonably close 
to 1 . The coefficient yj irv/2 in the Eq. (J7J) is partly con- 
nected with normalization of the partial radial functions 

ifii(kr) 



2 V 2 '- 1 [Mm 2 , 



The ratio of the reduced widths calculated at square well 
boundary by Eq. ([5]) and for " effective infinity" by Eq. 
([5]) is approximately 

el/el « (21 - i)/(2i + 1) , 

again the expression becomes an identity for E r — * 0. 

In Table Q] the mentioned values are provided for the 
simple model employed here. We can find out that (i) 
the effect of the subbarrier region is very important for 
weak barriers (e.g. I = 1 in our example) and it gradu- 
ally diminish as the barrier grows; (ii) the effect is not 
absolutely negligible even for quite high barriers; (iii) the 
direct numerical comparison shows that the DRW value 
6*2 calculated for "effective infinity" is exactly consistent 
with the width obtained in scattering calculations, and 
should be considered as a correct definition. 



ipf(k'r)cpi(kr)dr = -S(k' - k) 



and partly with conversion from integration over dk to 
integration over dE. Now substituting (J7J into |T]) and 
using the identity 



T(E r ) = T(E r ) -a 



f*dri>f(k a ,r)(V-V) Mkr,r) 
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[E r -E a ) 2 +T a (E r f/4 



(8) 

is obtained. 

If the difference between energies E r — E a is small com- 
pared to the height of the barrier, the quasibound WFs 
ipl and ipi should be quite close to each other and provide 
close DRW values 

E r - E a « E bar , $t(ka,r) » Mkr,r) -» 6 2 a /0 2 « 1 . 

(9) 

Also the variation of the kinetic energy should be small 
in such a case and the whole variation in energy should 
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be due to the potential energy change and it should be 
true that 



e 2 



dr &(k a ,r)(V-V) i>i{k r ,r) 



{E r — E a ) 2 Bi 



(10) 

with coefficient Bj sj 1. The examples of actual calcu- 
lated values of these coefficients are provided in Table 
U for the model of the previous Section. They indicate 
that the auxiliary resonance E a should be sufficiently 
narrow and sufficiently close to the physical resonance 
E r to make the approximation of Eq. ([7]) really precise. 
In the opposite case the value D[ needs to be correspond- 
ingly renormalised to preserve the identity ([H]). Anyhow 
we keep in mind that for states broader than 10~ 2 MeV 
the substitution ([7]) provides results which are valid only 
within a factor of 2. 

The issues discussed above are not of importance for 
the standard R-matrix phenomenology, but they should 
be clearly understood before we turn to the three-body 



III. THREE-BODY CASE 

In paper [ll[ we used the simplified three-body Hamil- 
tonian in which the proton-proton FSI was neglected [28| . 
In this Hamiltonian we used also a simplified expression 
for the three-body Coulomb interaction, which allows an 
isolation of the degrees of freedom. The three-body width 
in this simplified model was obtained as 



r(E 3r ) 



8E 3r , M x M y , „ / N|2 



de 
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I" Jo k x (s)ky(e) 
where energy and momenta of the subsystems are 



E x = eE 3r , E y = (l- e)E 3r , h = y / 2M z E z 
The "energy distribution" coefficient is defined 



Me) 



dXdY V iJk x (e)X) cp ly (k y (e)Y) 



x V 3 (X,Y) VLU y s(kr,X,Y) 



(12) 



where the function (fLi^i si^n X, Y) is the radial part of 
the solution with the outgoing asymptotic behaviour for 
a simplified three-body Hamiltonian 



H - E 3r j (p=[H x + H y + V 3 {p) - E 3r j £ 
= [T x +T V + V X (X) + V y {Y) + V 3 {p) - E 3r ) (p = 



at the complex pole energy E 3r — E 3r + iT/2. The three- 
body potential V 3 (p), which is depending only on the hy- 
perradius p, is used to form the three-body resonance 
and control the resonant energy. The functions ipi^ , <pi 
are the continuum eigenfunctions of the separable auxil- 
iary Hamiltonian H = H x + H y and they are normalized 



by the asymptotic condition ([2]) where phase shifts are 
defined for the subhamiltonians 

(T x + V x (X)-E x )<p lx (k x X) = 0, 
(T y + V y (X)-E y ) V > ly (k y Y) = 0. 

In the case when there are narrow resonant states in 
both "X" (at energy E xa ) and "Y" (at energy E ya ) sub- 
systems, the substitution ([7]) can be used for both <p/ x 
and ifi y . It is clear that this approximation is physical 
only when the so-called "Y" Jacobi system is chosen. In 
such a Jacobi system X is a coordinate between the core 
and a proton, Y is a distance between the "X" subsystem 
center of mass and the second proton. The hyperradius 
is then defined as 

n A r „o A r -\- 1 



1 



X' 



Y z 



where A c is the core mass. This approximation is good 
enough only in relatively heavy nuclei where the Y Ja- 
cobi coordinate is close to the single-particle proton co- 
ordinate due to a large core mass. The total auxiliary 
Hamiltonian H has a resonant energy 



E a — E xa + E ya , k ia — \j2MiE ia . 
Denoting the following integral as 

POO 

(V 3 ) = / dXdY0l(k xa ,X)0t y (kya,Y)V 3 (p) 
Jo 

x ip L i m i y s(kr,X,Y) , 

where ip are normalized quasibound WFs for the subsys- 
tems, the width is obtained as 



(11) T(E 3r ) = 
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(E x 



Dy T ya (Ey) 



{Ey-E ya y+T va {E y y/A 



E xa ) 2 +T xa (E x )y4 
(13) 



The expression for the two-proton width was, for the first 
time, obtained in a similar form in Ref. [121] . It is difficult 
to say, why this work did not attract attention and why 
its results have never been used. A possible reason could 
be the relatively complicated procedure used in Ref. [l2j 
and a lack of qualitative investigation of the model prop- 
erties and the nature of approximations involved. 
Using the following notations 



Ti(Ei) =2«ftPl t (E i ,r ch i,Z i ) = 
V u Je,E) = P lx (eE,r chx ,Z x 



A h (e 7 E,E ia ) 



(1 



Mir 2 chi 

Eia 



Ph(Ei,r c hi,Zi) , 

— e)E, Tcfry, Zy) , 



T ia (eE) 2 
4E 2 



we can rewrite Eq. (113|) in the "dimcnsionless" form for 
the penetration part of the expression: 



T(E 3r ) 



(V 3 ) 2 / deV u (e,E 3r ) 



D X Dy 2j 2 X 2jyE 3[ 

2t El a El a 
A lx (s,E 3r ,E xa )A ly (l-e,E 3r ,E ya ) . (14) 
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The expression for width is now explicitly factorized into 
"preexponent" (which has the dimension of energy) and 
dimensionless "exponential" part 



deVij y A^A^ ~ exp 



- \j M/E 3r 2Z COIe C 



where C is coefficient of the order of unity. 

Following the discussion of Eq. (|TUJ) we can write 



(V3) 2 = (E 3r — E xa — Eya) 2 D 3 



(15) 



where the dimensionless coefficient D 3 is presumed to be 
close to unity. Finally we get for the width: 

r(p > _ D x DyD 3 2 J 2 y E 3r (E 3r - E xa - E ya f f 1 
[ 3r) ' 2. M x M v r 2 chx rl hy El a El a J d£ 

xV Uy (e, E 3r )A Li (s, E 3ri E xa )A ly (1 - e, E 3r , E ya ) (16) 

The D 3 values calculated in the simplified three-body 
model are given in Table [TTJ It can be shown that the ra- 
tio (V 2 ) I (V3) 2 is a measure of the WF to be outside the 
interaction region (this statement is trivial to check for 
the square well potential). The values of (V^ 2 ) and (V3) 2 
calculated in the three-body model for 17 Ne and 45 Fe are 
also provided in Table HT1 The ratios (V^ 2 ) / (V 3 ) 2 are rea- 
sonably consistent with the D 3 values and are quite close 
to unity. This indicates that the WFs are predominantly 
localized in the "internal" regions. 



IV. DISCUSSION 

Table IIIII demonstrates a sensitivity of three-body 
widths estimated by quasiclassical expressions to differ- 
ent ingredients of the models. However, before discussing 
these effects we should make some overview of the mod- 
els. 



for the preexponent choice was like follows. Let's con- 
sider the two-body case in the situation of no barrier (no 
Coulomb interaction and zero angular momentum 1 = 0): 



T = 2 1 2 Pi(k 1 r ch ,Z) 



2 (kr ch ) e 2 ^2 2v 
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r c h r{r ch /2) 



The width in that case is formally just the inverse flight 
time for distance r ch /2 [denoted as r{r ch /2) above]. Us- 
ing the same assumption (no Coulomb and zero angular 
momentuma l x — l y — 0) we obtain from Eq. (|17[) 



r s (E 3r ) 



7T(r chx r ch yf/^ My M x )^ J 



Using the integral value 
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the estimate for the width is obtained as 



r s (E 3 



3E. 



1/2 

3r 



2(r c/la; r c ^) 1 /2(M y M :c ) 1 /4 ^r ch T ^ ch ) 



Thus, the width in such a "no barrier" case is normal- 
ized to the inverse flight time to some "typical internal 
distance" r c h- 

The above idea of the preexponent derivation can be 
found too simplistic. However, calculations show that for 
the considered cases Eq. (|17[) demonstrates a good con- 
sistency with Eq. (|16p when the channel radii are chosen 
in a sound way (namely they should be close to the inner 
classical turning point of the Coulomb barrier [26|). 



B. Special cases of the present model 

Let's consider the some special cases of Eq. (TT 



A. Our previous quasiclassical model 



1. True three-body decay, E xa ^> Es r , E ya ^> E^ r 



In the pioneering work of Goldansky [1] the differential 
probability of the two-proton decay was estimated as 



w{E 3r ,e) ~ exp 



2-7T Z COTe i/M ( 1 



r E. 



3r 
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In our works @, H, H, 0, Hi the quasiclassical expression 
for the two-proton width was used in the form 



r s (E 3 



6E. 



1/2 

3r 



n{r chx r chy fl 2 (M y M x f/' 



■J deV llcly (e,E 3r ). 



(17) 

The "exponential" coefficient Vxj, (£7 E 3r ) is closely re- 
lated to the function w(E 3r ,e) above. The motivation 



In that case the main contribution to the energy in- 
tegral in (TT6")) is connected with e — 1/2. Replacing the 
slowly varying functions A by their constant values at 
£ = 1/2 we obtain 



Tm{E 3r ) 



E 3r (E 3r — E xa — E ya ) 



{E xa — E 3r /2) 2 (E ya — E 3r /2) 2 

- D x D y D 3l 2 xl 2 v [ deP Uv (e,E 3r ) . (18) 
n Jo 



Thus we obtain formula analogous to (|17[) . used in our 
previous works, but with different preexponent. It should 
be noted: 

(i) The preexponent in this form is explicitly depending 
on the resonance properties of the subsystems. 
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(ii) It can be seen in Table HTT1 that agreement is very good 
between the value T in Eq. (|16p and the approximation 
T m in Eq. (fl8|) . calculated neglecting the variation of 
functions Aj within the decay window. 

(iii) The dimensionless coefficient D x D y D 3 is not that 
different from unity (at the level of 10 - 50%). The D x 
and D y coefficients should be very close to unity for suffi- 
ciently narrow states, so we need to know only the prop- 
erties of the two-body subsystem to fix this ingredient 
of the model. However, the D 3 coefficient is beyond the 
R-matrix ideology and requires a validation within the 
three-body model. 

(iv) If both subsystems (on X and Y coordinates) are 
equivalent, then E xa = E ya and the three-body width 
dependence on the two-body resonance position should 
have the following typical dependence: 

T(E 3r )~(E xa -E 3r /2)- 2 . (19) 

This dependence we have observed in the calculations 
within the simplified three-body model for 45 Fe (two 
equivalent final state interactions), see Fig. 14 of Rcf. 



2. Sequential decay, E xa < E- ir , E ya » E 3r 

It is quite simple to integrate over de in the Eq. (fl3|) 
for these energy conditions: 

r(E 3r ) « — V— 3 ^ — -2 T ya (E 3r - E xa ) 

(£j 3r — tj xa — hj ya ) 

Y(E 3r ) ps D y D 3 T ya (E 3r — E xa ) 

Thus the width is reduced to a two-body expression with 
some modifications, which takes into account the three- 
body character of the model considered for the resonant 
state. 



It should, however, be kept in mind that the considered 
approximation is valid and precise if the energy of the 
higher resonance E ya is kept significantly below the bar- 
rier. It is necessary for two reasons, (i) Geometric prox- 
imity of the resonance WFs t/>j* ; (X, Y) and ijji x i v (X, Y) 
is requred for Eq. (fl~5|) which is based on the smallness 
of kinetic energy contribution to the variation of total 
energy, (ii) The width of the upper resonance should be 
sufficiently small as demonstrated in Section Ul Bl All the 
mentioned conditions for this special case are difficult to 
meet in practice. It can be seen in Fig. 5 of Ref. [ll| that 
the width calculated for the 17 Ne 3/2" state in OFSI 
model (one subsystem is nonresonant) follows the ana- 
lytical dependence Eq. (fT9|) with significant deviations, 
while in the case of calculations for 45 Fe shown in Fig. 13 
[liT | (both subsystems have resonances) the agreement is 
practically perfect. 

A quasiclassical expression for the width was intro- 
duced in the papers [H, HBl for 12 O in the form 

T b (E 3r ) = 2 _£±lk^L ffeV lxly {s,E)Aa£,E 3r ,E xa ) , 

7T E xa Jo 

(21) 

(given here in our notations), which is practically equiv- 
alent to Eq. (f2TJ)) . It was introduced without derivation 
and analysis of the involved approximations. The ex- 
pression was utilized in the series of works [l7], EH [H, 
|2(I I2H [2^ | by Barker and Brown in somewhat modified 
form. We can see in Table HTT1 that for the systems under 
consideration the disagreement between "correct" value 
r Eq. (fl6|) and "special case" r& Eq. (f2Tj) is as large as a 
factor 2 — 4. It is also easy to show analytically that for 
systems with E xa = E ya (that means for all the ground 
state decays) the difference should be about a factor 4. 



C. Diproton model 



3. True three-body decay, E xa > E% T , E ya 3> E$ r 

In this case we can approximate the (E 3r — E xa — E ya ) 2 
by E ya in the numerator and (E ya — E 3r /2) 2 by E ya in 
the denominator and obtain 



r(E 3r ) = - 



2 D x D y D 3l 2 xl 2 y E 3r 



it {E xa - E 3r /2f 



[ eV Uy (e,E 3r ). (20) 
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This situation is close to the one final state interaction 
(OFSI) model considered in Ref. [Ill for methodological 
purposes. The three-body width dependence on the two- 
body resonance position is again 

T(E 3r ) ~ {E xa — E 3r /2)~ 2 . 

as in Eq. (fT9|) where two equivalent FSIs are considered. 
This behaviour is reasonably well reproduced for 17 Ne 
within the simplified three-body OFSI model, see Fig. 5 
in Ref. [13 . 



The nature of approximations used for derivation of 
the quasiclassical formula (flB|) is such that they can not 
be used for derivation of a corresponding formula for a 
"diproton" model. The formula (|2"T1) [actu ally analogous 
formula] was used in the papers Hi M M M il m 
the "diproton" form by formally making the "X" subsys- 
tem to be subsystem of the p-p motion. We find serious 
problems in such an approach. 

We can consider this issue from a different side: the 
analysis can be started from the original expressions (|lip , 
([T2")) [not from approximation [)14p] and the situation 
studied with one resonant and one nonresonant subsys- 
tem. The direct integration in (Ti"2"|) for the case of one 
resonant and one nonresonant continuum in the subsys- 
tems leads in the general case to the expression analogous 
to pTjl . However it contains a complicated additional co- 
efficient that is strongly dependent on fine details of the 
system geometry and the behaviour of the nonresonant 
continuum. Thus we can conclude that in this case Eq. 
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(j2Tj) should be correct only within an order of the mag- 
nitude. 

The analysis provided in Ref. for the case physi- 
cally corresponding to the "diproton" model (the OFSI 
model in the "T" Jacobi system) showed that the situa- 
tion in that particular case is even worse. The semiana- 
lytical three-body model of Ref. [ll[ provides the width 
values which are very small (compared to those typically 
evaluated in diproton model). They could be matched to 
those obtained from Eq. (|2Tj) only if very small channel 
radii r c hy (1 — 2 fm) are chosen for emission of the "dipro- 
ton" . This requirement is certainly not consistent with 
the practice of the "diproton" model application, where 
the r c h y is typically chosen as a nucleus radius plus some 
"radius of diproton" . 



D. Relation to the three-body calculations 

Some important points could be understood using the 
information listed in Table HTT1 

Treatment of the Coulomb interaction in the present 
model follows the approximations discussed in detail in 
Ref. [llj: "no p-p" means that Coulomb interaction be- 
tween protons is just neglected (product of charges are 
Z x = Z y = Z core both in X and Y subsystems) ; effective 
treatment of Coulomb "Eff." means that Coulomb inter- 
action in Y coordinate is formed by proton and "effective 
particle" core+proton (product of charg es are Z x — Z core 
in X and Z y = Z core + 1 in Y subsystem). Sensitivity 
to the choice of the Coulomb treatment is relatively high 
(factor 1.5 — 5, depending on particular model). 

The values of T derived from Eq. (fT6|) reasonably well 
agree with the corresponding results of the three-body 
calculations with the simplified Hamiltonian from Ref. 
PH . The disagreement can be reduced if (i) the widths 
of the subsystems (r^, T y ) are taken the same as in the 
three-body model, (ii) corrections for coefficient D3 (see 
Table HT|) are taken into account and (iii) radii of channels 
in Eq. (I16p are close to the inner classical turning points 
[26| for corresponding potentials in the three-body model 
(see Table Hp. 

The sensitivity of the results to the "unphysical" (not 
leading to modification of observable values T xa , T ya ) 
variation of channel radii is moderate (about factor of 
1.5) except for the T s model Eq. (TIT)) from our previ- 
ous works. However, even this model is providing results 
consistent with Eq. (fl6|) if the channel radii are chosen 
close to classical turning point. 

From comparison of the simplified three-body model 
with the realistic three-body model in Ref. [ll| we can see 
that the calculations using the effective Coulomb interac- 
tion are reasonably close to realistic results and should be 
a preferable choice. However, even calculations with the 
effective Coulomb interaction differs from the realistic re- 
sults typically by a factor 1.3 — 3. Thus the calculations 
in the quasiclassical model presented here (which itself 
is an approximation to the three-body model with sim- 



plified Hamiltonian of Ref. [Tl|) are not a replacement 
for the realistic three-body calculations if a better than 
mentioned precision is requested. 



V. CONCLUSION. 

In this work we derived the quasiclassical ("R-matrix 
type") formula for two-proton decay widths. The preex- 
ponent coefficients arc defined and evaluated numerically 
using the simplified three-body model. The following as- 
pects of the obtained results should be emphasized. 

(i) The derivation is based on the three-body model with 
a simplified Hamiltonian [ll| , which omits one FSI (p-p) 
and treats another in an approximate way (one of core-p 
interactions). The first approximation can be justified by 
the weakness of the p-p interaction compared to core-p 
interaction and should be good in heavy systems. The 
second approximation is connected with the finite mass of 
the core and also becomes well justified in heavy systems. 

(ii) The derivation of the quasiclassical formula requires 
existence of narrow states in both core-p subsystems. 
This condition is also well satisfied for the ground states 
of relatively heavy systems. 

(iii) The derived formula is basically the same as that ob- 
tained by Galitsky and Cheltsov in Ref. [T3| by a some- 
what different procedure. As far as we started from 
construction and validation of the simplified three-body 
model [HJ we were able to define a precision of approx- 
imations used for transition to the quasiclassical model 
and define ing redients of the model which remained un- 
defined in [12j. 

(iv) The most important dependencies of the three-body 
width (|16p are fixed by observable properties of the sub- 
systems (positions and widths of the lowest resonances in 
the subsystems). However, there is "unphysical" sensi- 
tivity to channel radii when the observables for the sub- 
systems are fixed. The quasiclassical formula provides 
a good agreement with the simplified three-body model 
when the radii of channels in the subsystems are cho- 
sen to be close and outside of the inner classical turning 
points of the barriers. 

jv) The formula (21]), used in papers [H GJ, El S3 EH, 
|22| , is a special case of Eq. (jTBf . It can be obtained by 
formally assuming the energy of one of the states in the 
subsystems to go to positive infinity. Such an assumption 
is unphysical when both valence protons populate states 
with the same quantum numbers. Thus formula (j2~lj) is 
valid within a factor which can be as large as 4. 
(vi) The derived expressions can not be used in the 
form of the "diproton" model without introducing a large 
(above an order of magnitude) uncertainty. 

Having in mind the origin and scale of the uncertainties 
introduced by reducing the realistic three-body model 
to a simplified three-body model (discussed in details in 
[HI) and by reducing the simplified three-body model to 
the quasiclassical three-body model (considered in this 
work) we now get a basis for the appropriate (within 
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the limits of its reliability) application of this model for 
estimates of the two-proton widths. 
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TABLE I: Properties of the test " 2 n + 2 n" system (the reduced mass is equal to the neutron mass) for square well potential 
R = 4 fm and E r =0.1 MeV. The dimensionless reduced widths calculated for well boundary (8 2 ) and for "infinity" in the 
sense of Eq. (|6]) (el)- Values Bi [Eq. (|10|l ] calculated for different energies of a resonance in auxiliary Hamiltonian. The width 
r obtained by two-body scattering calculations [it exactly coincides with F defined by Eq. © for reduced widths 8 2 ]. 



1 


el 


el 


ell el 


Bl{Ea = \.\E r ) 


Bi(E a — 2E r ) 


Bi{E a = hE r ) 


T (MeV) 


1 


2.044 


0.671 


0.329 


0.81 


0.75 


0.58 


3.46 x 10" : 


2 


2.013 


1.205 


0.598 


1.04 


1.06 


1.11 


5.58 x 10~ 


3 


2.007 


1.432 


0.714 


1.007 


1.011 


1.03 


2.14 x 10~ ! 


4 


2.002 


1.557 


0.777 


1.0022 


1.0046 


1.013 


3.54 x 10~ : 



TABLE II: Parameters for quasiclassical approximation calculated in a three-body (TFSI) model Ref. [ill ]. 





E 


l x 


E 


ly 


Eya 


E a 




(vi) 


(Vi)/(V 3 f 


D 3 


17 Ne 


0.344 





0.535 


2 


0.96 


1.495 


1.714 


2.426 


1.414 


1.338 


45 Fe 


1.154 


1 


1.480 


1 


1.480 


2.960 


3.285 


3.991 


1.215 


1.098 



TABLE III: Width sensitivity to the treatment of the charges in the subsystems and to channel radii for fixed properties of the 
subsystems "X" and "Y". The three-body widths F s Eq. (Tl7j), F b Eq. J2TJ, Y m Eq. JT8J, T Eq. |16J) with D s = 1, are given in 
10" 14 MeV units for 17 Ne and 10" 19 MeV units for 45 Fe. The two-proton decay energies E 3r are 0.344 MeV for 17 Ne and 1.154 
MeV for 45 Fe. The recent experimental data on 45 Fe ;25j] provide F 2p = 2.85±g;^ x 10" 19 MeV \T 1/2 (2p) = 1.6±g;| ms] for 
Ezr = 1.154(16) MeV and two-proton branching ratio Br(2p) — 0.57. The E xa , E ya values are given in Table HU The widths 
T xa and T ya are chosen to be the same as in the corresponding potential model of Ref. Tables I, II. 



Coulomb 


I, 


e 2 


r c hx (fm) T. 


za (keV) 


ly 


el 


Tchy (fm) 


T ya (keV) 


T s 


r 6 


r 

J- m 


r 


17 Ne No p-p 





0.986 


3.53° 


17.9 


2 


1.85 


3.53 a 


3.5 


2.63 


1.75 


3.70 


3.72 


No p-p 





0.667 


5.06 6 


17.9 


2 


1.145 


4.04 6 


3.5 


9.44 


2.04 


4.33 


4.35 


Eff. 





0.667 


5.06 6 


17.9 


2 


1.145 


4.12 


2.2 


1.86 


0.385 


0.845 


0.844 


Eff. 





0.667 


5.06 6 


17.9 


2 


0.162 


7.0 C 


2.2 


2.27 


0.510 


1.12 


1.10 


Eff. 





0.363 


8.0 C 


17.9 


2 


1.145 


4.12 


2.2 


5.84 


0.522 


1.15 


1.13 


Eff. 





0.363 


8.0 C 


17.9 


2 


0.162 


7.0 C 


2.2 


71.2 


0.691 


1.52 


1.50 


45 Fe No p-p 


1 


1.07 


4.72 6 


0.257 


1 


1.07 


4.72 6 


0.257 


10.5 


3.24 


12.7 


12.9 


No p-p 


1 


0.888 


4.94° 


0.257 


1 


0.888 


4.94 a 


0.257 


16.6 


3.36 


13.1 


13.3 


Eff. 


1 


1.07 


4.72 6 


0.257 


1 


1.03 


4.76 


0.15 


3.94 


1.14 


4.54 


4.62 


Eff. 


1 


0.888 


4.94° 


0.257 


1 


0.85 


4.98 a 


0.15 


6.24 


1.18 


4.70 


4.78 


Eff. 


1 


0.143 


7.5 C 


0.257 


1 


0.128 


7.56 


0.15 


667 


2.02 


8.02 


8.16 


"This is radius from 
the papers 0,H&S 


8] 


'stematics 


r ch = 1.4(A C orc - 


1- l) 1 / 3 used 


in 

















6 This is radius of the classical inner turning point for potential 
Ref. [ll, Tables I, II. 

c At this radius the nuclear potential becomes negligible in poten- 
tial model Ref. [H, Tables I, II. 



